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Abstract 

We present new upper and lower bounds to the primordial abundances of deuterium 
and helium-3 based on observational data from the solar system and the interstellar 
medium. Independent of any model for the primordial production of the elements we 
find (at the 95% C.L.): 1.5 x IQ-^ < {D/H)p < 10.0 x 10"^ and {^He/H)p < 2.6 x 10"^ 
When combined with the predictions of standard big bang nucleosynthesis, these con- 
straints lead to a 95% C.L. bound on the primordial abundance of deuterium: [D / H)^^^^^ = 
(3.5l^'g) X 10~^. Measurements of deuterium absorption in the spectra of high redshift 
QSOs will directly test this prediction. The implications of this prediction for the primor- 
dial abundances of helium-4 and lithium-7 are discussed, as well as those for the universal 
density of baryons. 
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Introduction 



Of the primordially produced light elements D, ^He, "^He, and ^Li, the predicted 
abundance of D is the most sensitive to the baryon density (technically the baryon to 
photon ratio at the time of nucleosynthesis) and thus an accurate determination of the 
primordial D abundance serves, via standard big bang nucleosynthesis (BBN), as a sensitive 
diagnostic of the baryon density of the universe. Conversely, the primordial abundance of 
D can be used as a probe of the BBN model - any D abundance corresponds to definite BBN 
predictions for the abundances of ^He, "^He, and ^Li which can be tested with observational 
data. Unfortunately, D is the most weakly bound of the stable nuclei and is therefore 
easily destroyed in astrophysical environments (Epstein, Lattimer and Schramm 1976). 
This makes the extraction of the primordial D abundance from observations of the present 
interstellar medium (ISM) and from the solar system a complicated and, necessarily, model- 
dependent task. Initial studies (Yang et al. 1984 (YTSSO); Dearborn, Schramm and 
Steigman 1986 (DSS); Walker et al. 1991 (WSSOK); Steigman and Tosi 1992 (ST92)) 
focused on the evolution of D plus ^He in an attempt to soften the dependence on chemical 
evolution. Recently Steigman and Tosi (1995(ST95)) improved on this earlier work by 
introducing a generic parameterization of D and ^He chemical evolution. Here we show 
that the ST95 analysis leads to a bound on the primordial abundances of D and ^He 
independent of the mechanism by which these elements are primordially produced. 

Our constraints, when combined with the BBN D and ^He predictions, lead to a 
prediction for the primordial abundance of D which best satisfies the observational data 
from the ISM and solar system and the BBN predictions: {D j H)^^^-!. = (S.Sl^'g) X 10~^ 
(unless otherwise noted, quoted errors throughout correspond to 95% C.L.). This predicted 
abundance is roughly an order of magnitude less than that inferred from the possible 
detection of D in a QSO absorption line system by Songaila et al. (1994) and Carswell et 
al. (1994) but is in agreement with a possible D detection (Tytler 1994) in a different QSO 
absorption line system. For this narrow range of best fit values for primordial D, BBN 
predicts a ^He mass fraction of Y^,^,^ = 0.247 ±0.004 and CLi/H)^^^^ = {2.9t*^-l) x 10"^°. 
The predicted "^He abundance does not agree with the primordial abundance as inferred 
from metal poor HII regions, possibly indicating an underestimate of the systematic errors 
in that data set. The predicted ^Li abundance is consistent with the data from metal-poor 
stars in our Galaxy's halo. The baryon density (relative to the critical density) consistent 
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with this D prediction is Q^h^^^ = (0.07lo'o2)5 which is roughly a factor of two larger 
than earlier BBN estimates. 

The Observational Bounds on Primordial D and ^He 

The determination of the primordial D abundance from present day observations of 
the ISM and solar system clearly requires knowledge of the chemical evolution history of 
the Galaxy - the D abundance decreases with time since D is destroyed in astrophysical 
environments. Using the argument that some of the destroyed D ends up as ^He and 
some ^He survives stellar processing, past studies have focused on an upper bound to 
the abundance of primordial D plus ^He in an attempt to reduce the chemical evolution 
dependence. If the post big bang production of ^He is neglected, the dependence on 
chemical evolution is characterized by a single parameter, (/g, the fraction of ^He that 
survives stellar processing. The amount of ^He which survives stellar processing depends 
on stellar models, the initial mass function and the history of star formation, and therefore 
(/g hides all of our ignorance about the chemical evolution history of the Galaxy. Recently, 
Steigman and Tosi (1995)(ST95) provided a generic analytic framework for the evolution 
of D and ^He. Starting from this analysis we show below that it is possible to bound 
simultaneously both the primordial abundance of D and that of ^He. Any model which 
predicts the primordial abundances of D and ^He is subject to this constraint. Comparing 
our new constraints to the predictions of BBN permit us to predict accurately the best fit 
primordial D abundance. 

In order to trace the evolution of the D abundance we follow ST95 and track the 
fraction of gas, /(t), which never undergoes stellar processing by the time t. Since any D 
cycled through stars is destroyed, / and the D survival factor, are the same: f2{t) = 
X2{t) / X2 p = f{t) (we use X- to indicate mass fractions and the subscript "P" denotes 
the primordial abundance). Again following ST95, the mass fraction of ^He at any time t 
is bounded by: X^[t) > X^pf[t) + (1 — f[t))g^[X^ p + 3X2p/2), where is the fraction 
of ^He which survives stellar processing. Here the first term on the right hand side is 
unprocessed ^He, the second term is the amount of ^He (both primordial and secondary 
from primordial D burning) which survives stellar processing, and the inequality accounts 
for possible stellar ^He production.^ We can combine these expressions describing D and 

^ Recently, Olive et al. (1994) and Tosi et al. (1994) have examined models including 
net stellar production of ^He. 
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^He evolution to write a single inequality relating the the primordial abundances of D and 
^He {y2p and 1/3 p, the number fractions relative to hydrogen, respectively) to their solar 
system abundances (1/2 and j/g, respectively): 

1 / Vsp 



. X 

^ ^ X, 



9i \y23p 



^23 



where X is the hydrogen mass fraction [X/Xp = 0.92 it 0.08 (ST95)) and subscript "23" 
denotes D+^He. This is the equation discussed by ST95 and is an improved (due to 
the last term which corrects for primordial ^He) version of that used by YTSSO and 
subsequent analyses (DSS; WSSOK). The chemical evolution model dependence involved in 
this constraint lies entirely in the choice of (/g, the ^He survival fraction. ST95 examined the 
bounds to the baryon-to-photon ratio which result from eq.(l) if the BBN predictions for 
the primordial components are adopted. In order to constrain simultaneously primordial 
D and ^He, we may write eq.(l) as a quadratic in 1/2 p- 

X IX X / 1 \ 

y2p + (y3p - ^^2 — ^^3)^2 p + ^ — 1 y^y^ p-^- 
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Provided that y^p is not too large, there are two real solutions corresponding to the equality 
of eq.(2) which represent upper and lower bounds to 1/2 p consistent with the solar system 
abundances of D and ^He and some chemical evolution model for the stellar processing of 
^He. In addition we require that any particular allowed choice of y^p be greater than both 
the ISM and solar system abundances of deuterium (since deuterium is always destroyed 
in astrophysical environments, 1/2 p — ^2© — y2lSM where jsm ~ (1.6 it 0.4) X 10~^ 
from the combined Copernicus and lUE (McCuUough 1992) and HST (Linsky et al. 1993) 
observations of the interstellar medium D/H ratio and 1/20 = (2-6 it 1.8) X 10~^ is the 
solar system abundance of D(ST95)). In Fig.(l) we show the solution to eq.(2) in the 
y2p - y-^p plane using Geiss' (1993) data for the solar system abundances of D and ^He 
and Y = 0.28 it 0.04 for the solar "^He mass fraction (see ST95 and Thomas et al. (1995) 
for details). Following YTSSO and WSSOK, we adopt > 0.25 and show the 95% C.L. 
bound for (/g = 0.25 as a solid curve, the interior of which represents primordial abundances 
of D and ^He consistent with the solar system abundances of D and ^He. The contour is 
obtained from a joint likelihood function which includes gaussian probability distributions 
for all observed quantities and the assumption that for any particular choice of allowed 
1/2 p all values of j/gp are equally likely (that is, each point interior to our contour can 



be associated with some non-zero stellar ^He production, each value of which we consider 
equally likely since the ST95 parameterization carries no information regarding stellar ^He 
production)^. We also incorporate a rigorous statistical definition in treating upper/lower 
bounds [for a complete discussion of the primordial abundance error analysis, see Thomas 
et al. (1995) and Hata et al. (1995)]. 

For comparison, YTSSO, DSS, and WSSOK adopted (/g > 0.25 and ST92's results (as 
reanalyzed by ST95) show (/g 0.6. We feel our choice for (/g is conservative with respect 
to most chemical evolution estimates for (/g in the sense that recent studies find larger (/g 
than assumed here (see, for example, Olive et al. (1994)) and the fact for any star (/g > 0.25 
(DSS). The trend is that the larger (/g, the more restrictive our constraint. For example, 
if we choose (/g = 0.6 the bound on D for a given ^He abundance decreases by ~ 40%. We 
could also take (/g to be distributed between 0.25 and 0.5 with equal probablility. In this 
case, the upper bound to D for a given ^He abundance decreases by ~ 30%. In contrast, 
(/g could be smaller than 0.25 if a significant fraction of material has been cycled through 
two or more generations of stars. We have also plotted the bound from WSSOK on the 
primordial abundance of D+^He: 1/23 p — -'-^•S X 10~^ (this is the 2a bound recalculated 
with the new solar system abundances of Geiss(1993)). Our upper bound to primordial D 
is more restrictive than WSSOK due to our more rigorous treatment of the upper and lower 
bounds. Independent of any model for primordial D and ^He production and assuming 
chemical evolution consistent with > 0.25 we find, 

1.5 X 10"^ < (^^^ < 10.0 X 10"^ and (^^) < 2-6 x 10"^ (3) 

The Best Fit D and ^He Abundances 

The key to determining the best fit value of the primordial D abundance is to compare 
the BBN predicted D vs. ^He relation with our D vs. ^He bounds from solar system and 

^ We thank M. Turner for two useful criticisms of the original version of this paper: (1) 
we included only the ISM lower bound to y^p and (2) our original normalization scheme 
for the likelihood function assumed that each value of j/gp was equally likely and thus 
over- weighted the larger values of J/gp. The first correction increases the lower bound to 
D/H while the second correction only slightly alters the allowed region. 
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ISM data. In Fig.(l), we show the resuhs of a Monte Carlo simulation of the standard^ 
BBN reaction network. The dotted lines show the Icr bounds to the BBN predictions 
for D and ^He. The numbers interior to this band (3, 4, 5, ...) represent the values of 
riiQ{= f] X 10"^°, where r] is the baryon-to-photon ratio) corresponding to abundances of D 
and ^He predicted by standard BBN. Over the range of D and ^He abundances of interest 
the correlation of the errors in each abundance do not appreciably affect the BBN D vs. 
^He relation and we construct the la contour with the usual analysis. Our BBN code 
is an updated version of the Wagoner code (Wagoner 1973) with the errors in the reaction 
rates as detailed in Smith, Kawano, and Malaney (1993) (the exception being the neutron 
lifetime which we take to be = 887 it 2 s in accordance with the latest PDG (1994)) 
and the various corrections to Y^^^ as given by Kernan (1993). For further details, see 
Thomas et al. (1995) and Hata et al. (1995). Since D is burned more easily to ^He and 
than and ^He are burned to "^He, ^He decreases less rapidly with -q than does D and 
the y2BBN^^- ysBBN relation is monotonically increasing as seen in Fig.(l). Recall that 
high values of y2BBN Us BBN correspond to low values of rj and, vice- versa, so that in 
Fig.(l) r] increases from upper right towards lower left. 

As seen in Fig.(l), the constraint on D and ^He can be significantly narrowed when 
our model-independent constraint is combined with the BBN prediction. To determine the 
best fit abundances from a joint fit, we construct a likelihood function including both the 
observational distributions and the probability distributions of the BBN predictions. The 
combined fit region (shown shaded in Fig.(l)) is obtained from the A^^ method (PDG 
1994), where = — 21og£ and C is our likelihood function."^ This contour represents the 
best fit of the standard BBN predictions to our observational constraints. It corresponds 
to: ^ 

= (3.5l?:^) X 10"^ and (^) = (1.2 ± 0.3) x 10"^ (4) 

best V / best 

Discussion 

A clear confrontation between these BBN D predictions and anticipated observations 
of D in the spectra of QSO absorbers lies on the horizon. Recently Carswell et al. (1994) 

^ By "standard" we mean that the universe is assumed homogeneous and isotropic and 
that there are 3 light neutrino families. 

The method is justified when the likelihood function is a gaussian. Our likelihood 
function is well approximated by a gaussian in (y2P7 J/sp? v)- 
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and Songaila et al. (1994) reported {D I H)QgQ ~ (1.9 — 2.5) X 10""^ for a particular line of 
sight. This is roughly an order of magnitude larger than our prediction and at face value 
appears problematic for BBN and/or our treatment of chemical evolution. However, single 
QSO absorption system measurements may only provide upper limits to the D/H ratio due 
to the possibility of hydrogen interlopers masquerading as deuterium absorbers (Carswell 
et al. 1994; Songaila et al. 1994; Steigman 1994). In addition, reconciling such high D 
abundances with the ISM deuterium abundance and solar system ^He abundance proves 
difficult (Steigman 1994; Vangioni-Flam and Casse 1994; Olive et al. 1994). Recently, 
Tytler (1994) has reported a possible detection of D in a high redshift absorber at a level 
of ~ 2 X 10~^. While this apparently excellent agreement with our prediction could be 
viewed as a dramatic confirmation of standard BBN, it may indicate the natural dispersion 
of the inferred D abundances due to both hydrogen interlopers masquerading as deuterium 
and actual deuterium absorption at a level reduced by astration. Although it is clear that 
the QSO approach will provide an opportunity for direct measurement of the primordial 
D abundance, we caution the reader against drawing strong conclusions from the current 
statistically limited sample. 

There are indirect tests of the consistency our best fit deuterium abundance. For our 
D constraint, the corresponding best fit value of rj^^g^ = (S.Ol^'g) X 10""^°, is larger (by 
roughly a factor of two) than that which results from previous analyses(see, for example, 
WSSOK) which used ^Li along with D+^He to delineate the acceptable range of rj. We 
now examine the predictions of the other primordial light elements consistent with our 
best fit D abundance. From the data of Fig.(l), the predicted abundances corresponding 
to {DIH\^^, = {3.5tll) X 10-5 are: 

= (1.2±0.3) X 10-^ (5) 

best 

n,,, = 0.247 ± 0.004, (6) 

and 

(^)^ ^ = (2.91-^) X 10--. (7) 

The ^He prediction is currently the least testable of the three predictions because the 
chemical evolution of ^He is so poorly understood. This prediction is consistent with the 
smallest abundances of ^He inferred from observations of Galactic HII regions (Balser et 
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al. 1994) and supports the idea that the abundance of ^He increases as the Galaxy evolves 
even if net stellar production is neglected (ST95). 

The ^Li prediction is in good agreement with the primordial abundance of ^Li as in- 
ferred from abundances measured (see for example, Spite and Spite (1993) or Thorburn 
(1994)) in Pop. II halo stars : C LijH) 1.2 x 10""^° is the weighted average of the 
Li abundances for stars in the "Spite plateau" (WSSOK). The extraction of the primor- 
dial ^Li abundance from this data is model-dependent as the competition between stellar 
surface depletion and creation by cosmic ray nucleosynthesis must be taken into account. 
The standard depletion models (Deliyannis, Demarque and Kawaler 1990 and references 
therein) and the recent work of Vauclair and Charbonnel (1994) are consistent with modest 
depletion and an inferred primordial ^Li abundance of {J Lij H)p = (1.21q'5) X 10""^° (the 
errors take into account possible production and depletion (see for example Copi, Schramm 
and Turner (1995)), in good agreement with our best fit value (see Thomas et al. (1995) 
for further discussion). 

The largest discrepancy between the best fit predictions and observations occurs for 
"^He. This is not a new problem for BBN , but it becomes worse when one uses the predicted 
"^He abundance consistent with our best fit primordial deuterium result. Like all the other 
primordial elements, the primordial abundance of "^He must be inferred from data contam- 
inated by stellar processing. Linear extrapolation of the "^He abundance in metal poor HII 
regions to zero metallicity yields (Olive and Steigman 1995): Yp = 0.232 it 0.003 it /\Y^y^ 
where 0.003 is the Icr statistical error and AY^^^ represents all the systematic uncertain- 
ties in converting "^He"*" line strengths into "^He mass fractions. Estimates for AY^y^ range 
from 0.005 (WSSOK; Pagel 1993; Olive and Steigman 1994), to 0.015 (Copi, Schramm 
and Turner (1995)). Our best fit prediction for the abundance of "^He disagrees with the 
observational data unless the systematic errors are large 0.015j. If the systematic er- 
rors should not prove this large, then our best fit deuterium abundance challenges the 
consistency of the standard BBN model. ^ Such a potential conflict raises an interesting 
question: what is the best way to test standard BBN? We believe the answer crucially de- 
pends upon which of the primordial light elements can be best measured, both in terms of 

^ As theorists we would be remiss in not pointing out several variations to the standard 
BBN model which can resolve this crisis by reducing the predicted abundance of "^He: (1) 
massive tau neutrinos (Kawasaki et al. (1993); Dodelson et al. (1993)), (2) degenerate 
neutrinos (Beaudet and Yahil (1977)), (3) short-lived decaying particles (Scherrer and 
Turner (1988)). 
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the quality of the observational data and the model dependence of the extrapolation from 
non-primordial to primordial abundances. From our viewpoint, there are several choices: 
(1) use the QSO absorption system technique (or perhaps D absorption in the halo of our 
Galaxy and neighboring galaxies) to establish a nearly primordial D abundance, (2) reduce 
the systematic errors on measurements of Y in metal-poor HII regions, or (3) refine stellar 
models for ^Li depletion in Pop II stars. In option (1), the model dependence is small 
but the scatter in the data may be large. Both options (2) and (3) do not have the same 
baryon-diagnostic power (thru BBN) as does (1) and therefore the payoff from refining the 
data may not be as great. In two upcoming papers (Thomas et al. (1995) and Hata et 
al. (1995)) we discuss this in greater detail. 

Putting aside for the moment the possible conflict between our predicted and the 
observed "^He, the best flt D abundance corresponds to a baryon to photon ratio of f^i^g-,. = 
(5.0l^'5) X 10""^°. The upper bound to rj is model independent in the sense that it relies 
only on measuring the ISM and pre-solar abundances of D. The lower bound is more model 
dependent since it is based on an extrapolation of the primordial abundances of D and 
^He from solar system data. This range, higher by a factor of ~ 2 than that of WSSOK 
(due to our more restrictive constraint), corresponds to a baryon density parameter (for 
T^o = 2.726 ± 0.010 (Mather et al. 1990)): 

n^h,,' = (0.071°:°^). (8) 

For Hq < 100 km/s/Mpc, > 0.013 which reinforces the conclusion (WSSOK) that 
some, perhaps most, of the baryons in the Universe are "dark". Despite our higher range 
for r]^ the evidence for non-baryonic dark matter persists since for Hq > 40 km/s/Mpc, 
< 0.18. This higher range for rj goes part of the way towards resolving the x-ray cluster 
baryon crisis (White et al. 1993; Steigman and Felten 1994). Steigman and Felten (1994) 

1/2 

flnd that the current x-ray data suggest O/zgQ ' < 0.0731]^^ which for our best flt r] yields 

1/2 

the bound O/zgQ ' < 0.58 and so the x-ray "crisis" still exists in the sense that for = 1, 
Hq < 17 km/s/Mpc. 
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Figure Caption 



Figure 1. Shown are observational bounds, BBN predictions, and the best fit values 
for the primordial abundances of D and ^He. The 95% C.L. bound (solid curve) on 
the primordial abundances of D and ^He (by number relative to hydrogen) assumes a 
^He survival fraction (/g = 0.25. Primordial abundances of D and ^He interior to this 
curve are consistent with the solar system abundances of D and ^He and the assumed 
chemical evolution parameterization. Also shown are the interstellar medium bound (dash- 
dot curve) and the bound from earlier work (WSSOK) (dashed curve). The Icr range 
(dotted curves) of the BBN predictions for D vs. ^He (by number relative to hydrogen) 
are taken from our Monte Carlo of the BBN reaction network. The numbers interior to this 
band (3,4,5, ...) represent the values of rj^Q (see text for definition) corresponding to the 
primordial D and ^He abundances predicted by standard BBN. The shaded region (dotted 
curve) represents the 68 (95)% C.L. region for the abundances of D and ^He obtained by a 
comparison of the constraints derived from solar system and ISM observations (solid curve 
and dash-dot curve) and the predictions of BBN (dotted band). 



11 
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